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O . ABSTRACT 

^ ' We study the correlation between the 7-ray flux (i^-y), averaged over the first 1 1 months 

of Fermi survey and integrated above 100 MeV, and the radio flux density (F-^ at 20 GHz) 

Q>^ ■ of Fermi sources associated with a radio counterpart in the AT20G survey. Considering 

the blazars detected in both bands, the correlation is highly significant and has the form 
F^cxFr'^-*^*°°'*, similar for BL Lac and FSRQ sources. However, only a small fraction 

|jy ' (~1/15) of the AT20G radio sources with flat radio spectrum, are detected by Fermi. To un- 

derstand if this correlation is real, we examine the selection effects introduced by the flux 
limits of both the radio and 7-ray surveys, and the importance of variability of the 7-ray flux. 
After accounting for these effects, we find that the radio/7-ray flux correlation is real, but 



X 

II , its slope is steeper than the observed one, i.e. F-yCcF^'' with 5 in the range 1.25 — 1.5. The 

O ' observed F^-F^ correlation and the fraction of radio sources detected by Fermi is reproduced 



assuming a long term 7-ray flux variability following a log-normal probability distribution 



c/3 . with standard deviation a ^0.5 (corresponding to F^ varying by at least a factor 3). Such a 

^1 ' variability is compatible, even if not necessarily equal, with what observed when comparing, 

for the sources in common, the EGRET and the Fermi 7-ray fluxes (even if the Fermi fluxes 
PsJ ■ are averaged over '--^l year). Another indication of variability is the non detection of 12 out 

^ ' of 66 EGRET blazars by Fermi, despite its higher sensitivity. We also study the strong linear 

correlation between the 7-ray and the radio luminosity of the 144 PillQG-Fermi associations 
J;0 i with known redshift and show, through partial correlation analysis, that it is statistically ro- 

bust. Two possible implications of these correlations are discussed: the contribution of blazars 
to the extragalactic 7-ray background and the prediction of blazars that might undergo ex- 
t^^ ■ tremely high states of 7-ray emission in the next years. 

o 

(•^ I Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 

thermal — gamma-rays: theory — X-rays: general — radio continuum: general 
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1 INTRODUCTION Lacs (54), FSRQs (112), candidate blazars of unknown class (46) 

or other type of AGNs (10). Among the generic class of "AGN" 



The Large Area Telescope (LAT) on board the Fermi satellite (At- 
wood et al. 2009) detected 1451 sources (IFGL catalogue) in the 



there are different source types: starburst galaxies (NGC 253, M82 

- Abdo et al., 2010b), starburst/Seyfert 2 (NGC 4945 - Lenc & 
7-ray band above 100 MeV with a significance ^4.5(7 during Its „• '^nnn\ t t- ■. .■ i-nr j- 1 ^nx^o n^oc tc 

' ■', , ,, , , ^„,„ . ,„ , r ^ n^, Tingay, 2009) a Low-Excitation FRI radio galaxy (PKS 0625-35 

first 11 months survey (Abdo et al. 2010, AlO hereatter): 831 out „,. . ^ , ^„„„, .,■ it- ■. .■ t-tit j- 1 ^^ a 

m - Gliozzi et al., 2008) a High-Excitation FRI radio galaxy (Cen A 

of 1451 are classified as AGN (Abdo et al. 2010aU ^^^^ ^^ ^, ^^,^^ ^^,^„, ^ ^,^,._,, , .^^ ^^,^^_^ , ,p^^ ^^^^ 

We cross correlated (Ghirlanda et al. 2010, GIO hereafter) the 
Fermi IFGL catalogue with a complete flux limited sample of ra- 
dio sources detected by the Australia Telescope Compact Array 



Abdo et al. 2010f, 2010g), a Narrow-Line Seyfert 1 (PKS 2004- 

447 - Abdo et al., 2009b). The cross correlation of GIO also led to 

find 8 new associations among which two are classified as FSRQ 

and one as BL Lac. Therefore most of the 1FGL-AT20G associa- 
(ATCA) in a survey conducted at 20 GHz with a flux density limit . ,, , , r^cr,^ n, r 1 

„ „ tions are blazars of the FSRQ or BL Lac classes. 



The 230 1FGL-AT20G associations also have typically flat ra- 



^40 mjy (Murphy et al., 2010). The cross correlation led to identify 

highly probable (association probability >80%) radio counterparts 

f TTnicr^T ■ ,u 1C/-T AToar^ -^^ u ft dio spectra with spectral index in the range —0.5 < Q:(5_20 GHz) < 

for 230 IFGL sources, I.e. the 1FGL-AT20G associations hereafter ^ ^ r^ o "^ •' 

, 1 HI u » 1 oain\ ooo <■ »u i j i ■<; j ■ 0-5 and centred at Q(5_2oghz) ~ (with F^ cc ly ). However, 

(see also Mahony et al. 2010). 222 of these are already classified in ^° ^uv^nzj 

.1 j; . T AT- A /--XT ^^1 /iTA,- Auj ^1 iMir, N TIT tfic Tadlo AT20G soufces associated with a Fermi source of the 

the first LAT AGN Catalogue (ILAC - Abdo et al. 2010a) as BL 

IFGL catalog are only a minor fraction (--^1/15) of more than 3600 

AT20G sources with flat radio spectra (i.e. Q(5_2o ghz) > —0.5). 

* Email: giancarlo.ghirlanda@brera.inaf.it The 230 1FGL-AT20G associations show a statistically sig- 

^ http://heasarc.gsfc. nasa.govAV3Browse/fermi/fermilac. html nificant correlation between the 7-ray flux and the 20 GHz flux 
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density: F^oaFr^'^^ ^'^'^ (see GIO). This correlation has a similar 
slope when considering BL Lacs and FSRQs separately. 

The relevance of the F^~Fr correlation is twofold: it can help 
to estimate the contribution of blazars to the 7~-ray background (e.g. 
Stecker et al. 1993) and it could shed light on the physical link be- 
tween the emission processes in the radio and 7-ray energy bands. 
Indeed, the so called "blazar sequence" (Fossati et al. 1998, Ghis- 
ellini et al. 1998) was built by dividing blazars into bins of radio 
luminosity, thought to be a proxy for the bolometric one, and estab- 
lishes a link between the radio and the 7-ray emission. On the other 
hand, the radio and 7-ray emitting regions are probably different, 
since the rapid variability of the 7-ray flux suggests a compact size 
(see e.g. Tavecchio et al. 2010), for which the synchrotron spectrum 
is self-absorbed up to hundreds of GHz. We therefore believe that 
the link between the radio and the 7-ray emission (if any) must be 
indirect. One possibility is that both track the jet power, with the 
radio averaging it on a larger timescale than the 7-ray emission. 

The radio/7-ray properties of EGRET blazars suggested the 
possible existence of a Ky-Fr correlation considering bright 7-ray 
sources (Taylor et al. 2007) although at low radio fluxes no con- 
clusive claim could be made (e.g. Mucke et al. 1997). A possible 
correlation of Fj with the radio flux at 8.4 GHz was found in the 
population of blazars detected by Fermi in its first three months sur- 
vey (LBAS sample - Abdo et al. 2009). This correlation was more 
significant for BL Lacs (chance probability P = 0.05%) than for 
FSRQs (P = 8%). An updated version of the F^-Fi correlation, 
based on the Fermi first year AGN catalog, is reported by Giroletti 
et al. (2010). Recent studies of the radio-7 flux correlation in the 
LBAS sources (Kovalev et al. 2009a; 2009b) was conducted using 
the MOJAVE sample of extragalactic sources (with a flux limit of 
1.5 Jy at 15 GHz). Kovalev et al. (2009a) find that the parsec-scale 
radio emission and the 7-ray flux are strongly related in bright 7- 
ray objects, suggesting that Fermi selects the brightest objects from 
a flux-density limited sample of radio-loud sources. 

The F-f-Fr correlation is subject to the biases related to the 
flux limits of the radio and 7-ray surveys. However, these biases 
acting on the radio and the 7-ray surveys are independent. The 
AT20G-1LAC associations have been found by cross correlating 
two independent surveys: the radio AT20G (Murphy at el. 2009) 
and the Fermi 1 1 months survey catalogue (Abdo et al. 2010). 

The two main problems we want to tackle in this paper are: (1) 
understand why only a minor fraction (~1/15) of radio sources (of 
the AT20G survey) are detected in the 7-rays by Fermi, despite the 
possible existence of a correlation between the radio and the 7-ray 
flux; (2) recover the true radio-7-ray flux correlation by accounting 
for the selection effects of both the radio and the Fermi survey. 

An important aspect which could impact on these issues is the 
7-ray variability of blazars. We will consider in this paper two pos- 
sible variability patterns: a long term variability which is observed, 
for instance, when comparing the fluxes measured by EGRET and 
(almost ten years later) by Fermi for the sources in common, and 
a short term variability observed so far on daily timescales in the 
brightest sources. 

The paper is organized as follows: in §2 the correlation be- 
tween the 7-ray and the radio flux of the PJ'IOG-Fermi associa- 
tions is presented and in §3 the duty cycle of blazars in the 7-ray 
band is discussed. In §4 we describe the method used to reconstruct 
the true 7-ray radio flux correlation and in §5 we show our results 
and discuss their main possible implications. Summary and conclu- 
sions are given in §6. 




Log(F(20 GHz)) [mjy] 

Figure 1. The Fj—F^ correlation found with the IFGL sources with a 
counterpart in the 20 GHz radio survey AT20G catalogue (open squares). 
The solid line is the fitting correlation (with slope ~0.8). Also shown are the 
two main selection effects that we consider in this paper (see text): the radio 
flux limit (SI) at 40 mJy represented by the grey vertical shaded region and 
the Fermi detection sensitivity which increases for increasing 7-ray fluxes 
and it is represented by the bottom shaded region. The two hatched triangles 
represent the region of the plane where we should expect to find sources if 
there is no F^-F^ correlation. The absence of sources in these triangular re- 
gions is an indication of the existence of such a Fj-Fi congelation. The 3rd 
catalog EGRET sources (Hartmann et al. 1999) detected also by Fermi and 
not present in the AT20G survey (because at declination >0°) are shown 
with open (cyan) stars, those present in the AT20G survey by filled (green) 
squai'es and the 12 3EG sources not detected by Fermi are shown with (red) 
open circles. 



2 THE OBSERVED RADIO-GAMMA CORRELATION 

Fig. [T] shows the correlation found by GIO with the 230 IFGL 
Fermi sources associated with a radio counterpart in the AT20G 
sample (squares in Fig.[T). For these associations we computed the 
Fermi flux by integrating their spectrum (reported in the IFGL cat- 
alogue) above 100 MeV. The best fit (computed with the bisector 



method) of the Fj-Fr correlation is Fj(xFr 



. In the follow- 



ing we will refer to the latter as the "observed correlation" F^-F^. 
Our aim is to account for the possible selection effects acting on 
the F-j~Fi plane and recover the "real correlation" F-y-Fi which 
can have a different slope and normalization with respect to the ob- 
served one. With the "hat" quantities we indicate the 7 and radio 
flux of the sources generated through the simulations described in 
§4. To these sources we apply the selection effects in order to re- 
produce the observed correlation F^-Fi.ln the simulations (see §4) 
we consider the real radio sources adopting their real radio fluxes 
(reported in the AT20G survey). Therefore, F^ and _Fr coincide, 
whereas for each real radio source, the F^ is that obtained by as- 
suming a certain 7 ray variability (described in §3). 

The two main selection effects that could bias the F^-F^ cor- 
relation are the flux limit of the AT20GHz survey F^ >40 mJy (S 1) 
and the Fermi detection efficiency in the O.I-IOO GeV energy band 
(S2). These are schematically shown in Fig.[T] While the AT20G 
radio survey has a well defined 20 GHz flux limit (shaded grey re- 
gion in Fig.[TJ, the Fermi sensitivity depends on several parameters 
(Abdo et al. 2010c) like the source spectrum in the GeV band and 
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its position in the sky where the different intensity and anisotropy of 
the galactic and extra galactic 7-ray backgrounds can limit the de- 
tection efficiency as a function of the source flux. For instance, sim- 
ulations of sources distributed at high galactic latitudes \b\ > 20° 
show that the detection efficiency is only 1% for sources with F^ 
^ 10~* phot cm^^ s^^ although such a flux is above the lowest 
flux measured by Fermi (Abdo et al. 2010c) . This is due to the com- 
bination of two main effects: the intensity of the local background 
and the source spectral index (see Abdo et al. 2010c for details). We 
show in Fig. [T] the Fermi detection efficiency as a shaded coloured 
region. This is obtained from the efficiency curve reported in Fig. 7 
of Abdo et al. (2010c) that assumes a distribution of 7-ray photon 
spectral indices centered at P = 2.4 and with a dispersion of 0.28. 
Considering the distribution of sources in Fig. [T] and the two 
instrumental selection effects S 1 and S2, we note that there are two 
regions in the F^-Fi plane (the hatched regions in Fig. [T) where 
there are no sources. Although source number counts decrease with 
increasing fluxes (both in the radio and 7-ray band), there seems 
to be no apparent instrumental selection effect preventing the de- 
tection of bright 7-ray sources with intermediate/low radio fluxes 
(in the yellow-hatched upper left triangle) as well as bright radio 
sources above the Fermi detection limit (in the green-hatched lower 
right triangle). This suggests that the observed correlation is true 
although we expect that its real slope and normalization can be dif- 
ferent from those derived from the observed sources in the F-,-Fi 
plane because the latter is strongly biased at low F-, and Fr by the 
instrumental selections effects SI and S2 (as shown in Fig.[T). 



3 THE DUTY CYCLE OF THE 7-RAY FLUX IN 
BLAZARS 

The variability of the 7-ray emission in blazars has been since long 
discussed in the literature. It is now with Fermi that robust claims 
can be done, thanks to an almost continuous monitoring of the 7- 
ray sources in the sky with a relatively high sensitivity which allows 
to probe both the variations of the flux on long timescales (increas- 
ing with the mission elapsed time) and on short timescales (from 
months down to days for the brightest sources). 

Variability of the 7-ray flux of blazars could be the key in- 
gredient to explain why only a minor fraction of the radio sources 
detected in the AT20G survey have been detected by Fermi in the 
7-ray band and it could help to reconstruct the true F-,-Fr correla- 
tion accounting for the Fermi detection selection effect. 

In the following we will refer to two main variability patterns 
of the 7-ray flux of blazars: (1) a long term variability which seems 
to follow a log-normal probability distribution with a ~0.5 and (2) 
a short term variability which follows a non symmetric probability 
distribution skewed towards low flux levels. 



3.1 Long-tlmescale variability 

In the 3rd EGRET catalog (3EG - Hartmann et al. 1999) there are 
66 high confidence AGNs. 54 of these are detected by Fermi, while 
12 are not present in the ILAC catalog. Among the 54 BEG sources 
detected by Fermi, 48 have a published radio flux density at ~ 20 
GHz. Those in the northern emisphere can be added to Fig.[T](open 
cyan stars) and those in the southern emisphere (already present 
in the 1LAC-AT20G associations) are highlighted in Fig[T] (filled 
green squares). 

The 48 3EG sources (out of the 54 detected by Fermi), for 
which we could find the radio flux density (Fig. \Q, are consistent 



with the Fj-Fr correlation found through the 1FGL-AT20G asso- 
ciations. 

The 12 3EG sources (classified as blazars in Hartman et al. 
1999) not detected by Fermi in its 11 months survey are shown in 
Fig.[T]using their 3EG 7-ray flux (open red circles). 

The EGRET flux of these 12 sources is above the Fermi detec- 
tion sensitivity (shown by the shaded region in Fig.[Tll. They occupy 
a region where no apparent instrumental selection effect is present. 
Therefore, the non detection of these 3EG sources by Fermi must 
be due to their 7-ray variability over about a decade. 

It is interesting to compare the 7-ray flux of the sources de- 
tected both by EGRET and (about 10 years after) by Fermi. When 
doing this, we must recall that the Fermi fluxes are averages over 
the 1 1 months of the survey, while the EGRET fluxes corresponds 
to averages over a shorter time interval, typically few months, since 
they are derived from pointed observations. 

Fig.[2]shows the 3EG fluxes and those measured by Fermi for 
the common sources. We also show in Fig.|2]the 12 3EG sources not 
detected by Fermi as upper limits (green arrows in the top panel) 
and as lower limits on the EGRET to Fermi flux ratio (arrows in the 
bottom panel of Fig. |2j, obtained assuming, for Fermi, a limiting 
flux of 2.5x10^* ph cm^^ s^^, i.e. corresponding to the upper 
boundary of the shaded region S2 shown in Fig.[T] On average (see 
the bottom panel of Fig.|2j the EGRET fluxes were larger than those 
of Fermi. We find that the ratio of the flux measured by EGRET 
and by Fermi is distributed as a log-normal with standard deviation 
(7=0.5 (bottom panel of Fig.|2j. We note that this distribution fully 
comprises also the lower limits of the 12 3EG sources not detected 
by Fermi. However, given the presence of these lower limits in the 
bottom panel of Fig. |2] we also tested in our simulations a log- 
normal variability with (7=0.78 (represented in Fig.|2]by the dashed 
grey line), i.e. corresponding to a variation of the flux by a factor 
~6. 

If this is representative of a decadal flux variability of these 
sources, it explains why a fraction of EGRET sources were not 
detected by Fermi, despite its better sensitivity. 

This result, i.e. the possibility that the GeV flux of 7-ray 
blazars (even if time averaged over ~ 1 year) can vary by a factor 
3 (at la) over ~10 years will be used in the next section to recon- 
struct the true F-,-Fi correlation by accounting for the selection 
effect S2. 



3.2 Short-timescale variability 

Tavecchio et al. (2010) (see also Foschini et al. 2010) studied 
the variability of the GeV flux in the two blazars 3C 454.3 (see 
also Bonnoh et al. 2010) and PKS 1510-089 detected by Fermi 
and found variability on few days timescale by considering their 
emission as observed by Fermi in one year. During exceptionally 
bright events, significant variability was found also on intra-day 
timescales (Tavecchio et al. 2010; Abdo et al. 2010e). Tavecchio et 
al. (2010) also found that, in these two sources, the differential flux 
curve (representing the number of days a source spends at a given 
flux level F-j) has a similar pattern (also present in other sources 
- Tavecchio et al., in preparation) with a rising power law and 
a faster decay (steep power law) bracketing a characteristic peak. 
This short-timescale variability can be described as: 



N{F.y) (X 



(f7/-F7,brcak)° ' exp(-F-,/Fc, 
1 + (F^/Ky,brcak)''+" 



(1) 



where f-y. break is the flux corresponding to the break between the 
low-flux power law with slope a and the high-flux power law with 
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Figure 2. Top panel: comparison of Fermi and 3EG fluxes (> 100 MeV) of 
the 48 3EG sources detected by Fermi (reported in the ILAC sample) and 
for which we could find a radio flux density at ^^20 GHz in the literature. 
The 12 3EG sources not detected by Fermi are shown as upper limits (green 
arrows). FSRQs, BL Lacs and sources classified as "AGN" in the ILAC 
catalog are shown with different symbols (filled circles, open circles and 
open squares, respectively). The dashed line represents equality. Bottom 
panel: distribution of the ratio between the flux measured by EGRET and 
by Fermi for the sources in common. The 12 3EG sources not detected by 
Fermi are represented as lower limits (arrows). The hatched distribution can 
be represented by a Gaussian function (solid line) with central value /i=0. 17 
and standard deviation cr=0.5. It is also shown a Gaussian (dotted grey line) 
with a = 0.78, i.e. a factor 2 larger in linear scale, which is tested in the 
simulations (see §5). 



slope h and Fcut is the flux of the exponential cutoff. Since we 
are concerned with 7-ray fluxes averaged over one year, the short- 
timescale variability can produce very modest variations of the av- 
eraged flux. 

Instead, the long-timescale variability observed in the com- 
mon EGRET/Ferm; sources can change the flux by a factor 3 (at 
Icr level), implying a larger spread of the 7-ray flux, although cor- 
responding to longer timescales. 



4 THE SIMULATION 

We want to constrain the normalization and slope of the real F^- 
Fr correlation that reproduces the observed distribution of the 230 
1LAC-AT20G associations shown by squares in Fig. [T] and, at the 
same time, accounts for the non-detection of the large majority of 
the radio sources of the AT20G sample. Even if the 7-ray detection 
rate approaches 100% at the largest radio fluxes and decreases for 
lower radio fluxes, this is not a trivial task. This is because there are 
many radio sources, undetected by Fermi, with a radio flux compa- 
rable or even larger than those that are instead detected in 7-rays. 

By accounting for the Fermi detection efficiency and for the 
assumed duty cycle of blazars in the 7-ray band, we search for the 



F-y-Fr correlation that produces a distribution of simulated sources 
in the F^-F^ plane which matches the observed one. 

In particular, we consider all the radio sources with flat ra- 
dio spectrum in the AT20G survey and assign to them a 7-ray flux 
according to a given F-y-Fr correlation (whose normalisation and 
slope are the free parameters that we want to constrain). Then we 
shuffle the 7-ray flux of each source according to a law which is 
representative of the 7-ray variability and apply the Fermi detec- 
tion sensitivity to identify those simulated sources that should be 
detected by Fermi. In this way we populate the F^-Fi plane with 
simulated sources observable by Fermi for any assumed F-y-Fi cor- 
relation. We constrain the slope and normalisation of the F-,-Fi 
correlation by requiring that (i) the number of sources that should 
be detected by Fermi is consistent with the real number of sources 
defining the F-,-Fr correlation (i.e. 230±15), (ii) that the distribu- 
tion of Fi and F-, of the simulated sources is consistent with those 
of the real sources (this is evaluated through the Kolmogorov- 
Smirnov test, performed separately, between the Fi and F-, distri- 
butions of the simulated and of the real sources). The details of the 
simulation and its assumptions are described below, and we present 
our results in § 5. 

The simulation relies on some input assumptions: 



(i) the true F-y-Fi correlation is modelled as a power law 

log{F^) = K + 5log{Fr) 



(2) 



where the normalisation K and the slope 5 are the free parameters 
that we aim to constrain. The normalization is computed at 150 mjy 
in our simulations. This particular value corresponds to the average 
of the radio fluxes of the AT20G sources with flat radio spectrum 
used for the simulations; 

(ii) we consider the 3686 radio sources of the AT20G survey 
with flat radio spectrum, i.e. a{5_20GHz) > —0.5, similar to the 
radio spectrum of the 230 1FGL-AT20G associations defining the 
observed F-,-Fr correlation; 

(iii) the 7-ray flux variability: we assign to each radio source 
with a certain F^ (given by Eq.|2]l a flux F^ according to one of the 
two possible variability functions described in §3. First we assume 
the short-timescale variability function, described by Eq[T] Since 
we are concerned with averaged (over 1 1 months) fluxes, we extract 
11 7-ray fluxes from the N(F) distribution of Eq. [2] after having 
fixed its parameters to Fy. break ~ F^ and always setting a = 1.5, 
6 = 3 and Fcut = 5 Fy, break- We then average the 1 1 values of F-^ 
obtained in this way for each simulated source. The obtained flux is 
different from the initial F^, but by a small factor, and we anticipate 
that the dispersion induced by this treatment of variability is much 
smaller than the dispersion of the real sources in the F-,-Fi plane 
along the F^ axis. 

For this reason we adopted, as a second choice, the long-term 
variability function, i.e. a log-normal distribution with assigned 
standard deviation. This assumption is motivated by the compar- 
ison of the EGRET and Fermi flux for the common sources shown 
in Fig. [2] (see §3). The bottom panel of Fig. |2] shows that the dis- 
tribution of the flux ratio has a standard deviation of 0.5 which we 
also assume in our simulation. However, we will test a log-normal 
variability function with larger/lower standard deviation. 

We stress that while the variability function of Eq. [T] is repre- 
sentative of the blazar activity over one year, the log-normal with 
a — 0.5 found in §3 corresponds to a variability over ten years at 
least, i.e. the time between the EGRET and the Fermi measurement 
of the average flux of the sources in common. Furthermore, it is 
already indicative of how the ~1 year average 7-ray flux varies. 
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Figure 3. Example of a simulation of sources following the true correlation Fj-Fi (triple dot-dashed green line) with 5 = 1.3 and obtained with a 7-ray 
log-normal variability distribution with cr^ = 0.5. Filled black dots are the real AT20G sources with flat radio spectra (cj(5_20 GHz) > —0.5) simulated 
with the log-normal Fj variability function. The red points are 250 simulated sources which can be detected by Fermi according to its detection efficiency 
curve (selection effect S2 - Abdo et al. 2010c). The dot-dashed line is the best fit correlation of simulated sources and the solid blue line is the best fit to the 
real sources. The open blue squares are the 230 real sources of GIO. The bottom and right-hand side panels show the F^ and F~f distributions of the total 
simulated sample (dotted black line) and of the simulated sources after the application of the selection effects (red line). The distribution of the fluxes of the 
real sources detected by Fermi is shown by the blue line. 



not of the variations occurring on shorter timescales, as instead in- 
dicated by Eq.[T] For this reason, we will extract only one flux from 
this log-normal variability distribution. 

We did not model the possible radio variability of blazars. This 
is motivated by the fact that we simulate the 7-ray flux of real ra- 
dio sources, i.e. those with flat radio spectra in the AT20G complete 
survey. These are 3686 sources with Fr^ 40 mJy: their large num- 
ber ensures that we are sampling the possible range of variability 
of the radio flux density. 



4.1 One illustrative example 

In Fig. [3] we show an example of a simulation. For this example 
we have assumed a F^-Fi correlation with slope 5 = 1.3 and 
K = —8.6 at 150 mJy. After having assigned to each radio source 
its corresponding Fj, we have extracted its Fj from a log-normal 
distribution peaking at F^ and having a dispersion a = 0.5. The 
sources simulated in this way are shown by the grey dots in Fig. 
[3] and the assumed F^-F-r correlation is shown by the triple-dot- 
dashed (green) line. Then we selected the sources detectable by 
Fermi (red dots) with the following procedure. We have considered 



the Fermi 7-ray detection efficiency (selection effect S2) which is 
a function of the 7-ray flux (shaded area in Fig. [T). For S2 we 
used the detection efficiency curve presented in Fig. 7 of Abdo et 
al. (2010c) which was obtained through simulations of sources at 
galactic latitudes \h\ > 20°. 

For sources at low galactic latitudes the level of the galac- 
tic background and the larger number of sources may reduce the 
detection efficiency. Therefore, we mimic this effect by consider- 
ing a detection efficiency reduced by a factor of 3 for sources at 
—20 < b < 20°. This choice is motivated by Fig. |4] where it 
is shown the 7-ray flux distribution of the Fermi sources of the 
ILAC sample located at high galactic latitudes (|6| > 20° - solid 
histogram) and located along the Galactic plane (|fe| < 20° - dotted 
blue histogram). The two curves matching the left-hand side of the 
histograms represent the detection efficiency rescaled by a factor 3 
for the sources at low galactic latitudes. 

Among the simulated sources in a given bin of 7-ray flux we 
randomly extract a fraction of sources corresponding to the Fermi 
detection efficiency (from Abdo et al. 2010c) in that flux bin. This 
corresponds to the application of the S2 selection bias. The simu- 
lated sources surviving the S2 selection (i.e. "detectable sources", 
hereafter) are the red dots in Fig. [3] 
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Figure 4. Flux distribution oi Fermi sources of the ILAC sample located at 
high galactic latitudes (|fe| > 20° - solid histogram) and located along the 
Galactic plane (|fe| < 20° - dashed blue histogram). The two distributions 
are normalized to the corresponding solid angles. The solid and dotted red 
lines are two Gaussians fitting the flux distributions. The dashed (cyan) and 
solid (orange) lines show the detection efficiency (adapted from Abdo et al. 
2010c) and scaled by a factor 3 between the two histograms. 



Their distributions in the F^-F,. plane is compared with the 
distribution of real sources in the same plane. First, we compare 
independently the distributions of F^ and Fr of the detectable and 
real sources through the Kolmogorov-Smimov test and derive the 
corresponding probabilities P{KS)~, and P(KS)r- The two his- 
tograms (of the detectable and real sources) are shown in the two 
side-panels of Fig. [3] We consider that the detectable and the real 
sources have similar distributions in the F-y-Fr plane when the KS 
probabilities are both > 10^^. 

Then, for each assumed F^-F^ correlation (with fixed slope 
and normalization, 5 and K) we repeat the simulation 300 times 
and count the number of simulations yielding a KS probability 
>1% that both the radio and 7-ray flux distributions of detectable 
and real sources are drawn from the same parent population. 

We consider a set of input parameters (5, K) acceptable when 
more than 68% of the 300 simulations had P{KS).y and P{KS)r 
larger than 1%. 

Finally, among the acceptable simulations we identified those 
producing a number of detectable sources (red points in Fig. [3} 
equal to the real one (i.e. 230±15, open blue squares in the exam- 
ple of Fig. [3](. In the example shown in of Fig. [3] the number of 
simulated sources detectable by Fermi is --^250. These simulations 
give us the slope and normalization of the true F^-F^ correlation 
we are seeking. 



5 RESULTS 

In the following section we present the results obtained through 
our simulations under the two possible variability scenarios for the 
7-ray flux discussed in §. 3. 



5.1 Simulations with the short-term 7-ray variabihty 

In Fig. |5]we show the number of simulated sources detectable by 
Fermi versus the slope & of the F-y-F^ correlation. These are the 
results obtained assuming a short-term variability of the 7-ray flux 



described in §. 3.2. Each curve represents a different normaliza- 
tion K of the F^-Fr correlation. The open grey squares in Fig. \5\ 
are those simulations rejected because the distribution of F-y and/or 
Fr of the simulated sources are inconsistent with those of the real 
sources, i.e. P{KS)^ and/or P{KS)r < 10"^, in more than 68% 
of the repeated simulations. 

Vice versa, the filled circles correspond to distributions of sim- 
ulated sources (i.e. the red points in Fig. [3} in the F^-F- plane con- 
sistent with the distribution of the real sources, i.e. in more than 
68% of the repeated simulations (for each choice of the free pa- 
rameters 5 and K) the P{KS)^ and P{KS)r > 10"^ 

The results shown in Fig. [5] are obtained under the hypothesis 
that the short time variability of the 7-ray flux is described by Eq. 
[T] implying a very modest variation of the average flux. We note 
that all the acceptable simulations (filled circles) over predict the 
number of sources with respect to the 230 Fermi detected sources 
with an AT20G counterpart (GIO). Therefore a modest variability 
of the average 7-ray flux cannot reproduce the number of sources 



really observed in the F^- 
correlation. 



Fr plane for any assumed true F^-Fr 



5.2 Simulations with the long-term 7-ray variability 

The next step was to assume a larger variability function for F-y, i.e. 
a log-normal distribution with a^ — 0.5. The results of the sim- 
ulations under this assumption for the variability of F^ are shown 
in Fig. |6] (filled blue circles). In this case the acceptable F^-Fr 
correlations extend over a wider parameter range of normalization 
and slope, K and 5, of the F^-Fr correlation with respect to the 
acceptable simulations shown of Fig.|5] This is due to the assumed 
larger amplitude variability. In this case there is a set of simulations 
(those intersecting the horizontal shaded stripe in Fig. [6} that can 
also reproduce the real number of observed sources. Therefore the 
solutions we find in this case correspond to —8.7 < K < —8.3 
and 1.25 <S < 1.5. 

The results shown in Fig.|5]on the slope and normalisation can 
be understood with the aid of Fig.[3] The distribution of real sources 
in the F-y-Fr plane (open blue squares in Fig.[3l( is constraining. For 
instance, for very low normalisations K of the F-y-Fr correlation, 
the number of simulated points surviving the selection effect S2, 
i.e. the red points in Fig.[3l is lower than 230 and their distributions 
of F-y and F^ are inconsistent with that of the real sources. This is 
the case of the simulations below the shaded stripe (representing 
the number of real sources in the F-y-Fr plane) in Fig.|5] 

However, when a log-normal function is assumed for the vari- 
ability of F-y , the spread of the simulated points in F-y can be larger 
and, although the detectable sources are much fewer than the real 
one, their F-f and/or Fr distributions can still be consistent with 
those of the real sources. This explains why there are solutions 
in Fig. |6] which are acceptable although the number of detectable 
sources is smaller than the real one. 

Although the choice of a log-normal variability function with 
a-, — 0.5 is motivated by the long term variability of the EGRET 
sources detected by Fermi, we also verified how the solutions of 
the simulation depend on the choice of cr-y. In particular we tested 



a log-normal variability function with a^ = 0.18 and a^ 



0.78 



corresponding to a linear flux variation by a factor 1.5 and 6 re- 
spectively. The solutions are shown in Fig. [6] by the open star and 
circles, respectively. We find that u-y = 0.18 introduces a too small 
degree of variability (similar to the variability function of Eq. [TJ 
and all the solutions over predict the number of detectable sources 
with respect to the real number of associations. On the other hand 
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Figure 5. Simulation results. The number of simulated sources (which sur- 
vive the two instrumental selection effects SI and S2 described in the text) 
is plotted against the slope & of the assumed F^-F^ correlation. Each curve 
represents a set of simulations of the F-y~Fr correlation with fixed nor- 
malization K and varying slope 5. The open squares are those cases of (5, 
K) where more than 68% of the simulations fail to reproduce the observed 
source distribution of the real sources in the Fry-Fr plane (i.e. in more than 
68% of the simulations the simulated-source F-y and F^ distributions have 
KS probabilities < 1 % of being consistent with the Fj and F-y distributions 
of the real sources). The acceptable simulations are shown by the filled blue 
circles. The simulations are performed assuming Eq.[T]for the 7-ray vari- 
ability function that, after averaging, results in a very modest flux variabil- 
ity. The shaded region represents the number of real sources (i.e. 230±15) 
detected by Fermi with a radio counterpart which give rise to the observed 
F-y-Fr coiTelation. For reference, the vertical dashed line shows the slope 
of the observed F-i-F^ congelation. 
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Figure 6. Same as Fig.|5]but with a variability function for the 7-ray flux 
which is a log-normal with (7.^=0.5 (filled blue circles). Also shown are the 
results assuming a 7-ray log-normal variability function with (7.^=0.18 and 
a"-,=0.78 (shown by the open green stars and the open red circles respec- 
tively). 



a larger variability (i.e. a-y — 0.78) extends the space of acceptable 
solutions below those obtained with a.y = 0.5. Note that the po- 
sition of the red open circles and of the green open stars in Fig. |6] 
(corresponding to simulations with a-, — 0.78 and a^ — 0.15, re- 
spectively) does not exactly coincide with that of the filled blue cir- 
cles. This is because the three sets of simulations, shown in Fig. [6] 
have slightly different normalizations even when the slope is equal. 
For clarity, in Fig.|6]we draw only the curves (open grey connected 
squares) corresponding to the simulation with a-y = 0.5. 

As a caveat we stress that in our simulations, we adopted the 
Fermi sensitivity computed by Abdo et al. (2010c) which assume 
a spectral index distribution typical of FSRQ sources. While it is 
known that the Fermi sensitivity strongly depends on the source 
spectral index (e.g. Abdo et al. 2010a), it should be noted that 
our sample of 1FGL-AT20G associations is dominated by FSRQ 
sources. For BL Lac objects the better sensitivity of Fermi in de- 
tecting these hard sources would imply a lower detection limit (rep- 
resented by the shaded region S2 in FiglTJ. This would require, in 
order to reproduce the observed F^-F^ correlation, a combination 
of a slightly smaller normalization and slope of the F-y-F^ correla- 
tion possibly coupled with a slightly larger variability of the 7-ray 
flux. Still the results would be comparable with those derived with 
the detection sensitivity of FSRQ since we find that a variability of 
at least a factor 3 (i.e. (T=0.5) is necessary to reproduce the F-y-Fi 
correlation. 



5.3 Predicted number of FERMI detectable sources 

The Fermi sensitivity is increasing with the increasing exposure 
time of its survey. We can use the simulation and the reconstructed 
F-y-Fi correlation to infer the number of sources that will be de- 
tected with a future increase of the survey time which will improve 
the 1 1 months Fermi survey limit by a factor of 2. We expect that a 
larger number of detected 7-ray sources will have a counterpart in 
the AT20G radio survey. By running our simulation with the best 
correlation found in the previous section, we find that the total num- 
ber of southern sources present in ATG20G that will be detected 
by Fermi will go from the current 230 to ~430. This number is in 

agreement with what expected if the radio LogN—LogSr has slope 

" •^ 1.5 

-3/2 and considering the reconstructed correlation F-y oc Fr , and 

implies a 7-ray LogN — LogS-y with a slope flatter than euclidean. 



5.4 Gamma-Radio luminosity correlation 

The possible correlation between the radio and the 7-ray luminos- 
ity has been studied in the past with the aid of EGRET detected 
sources. Different groups reported a significant correlation between 
the radio luminosity (at frequencies larger than 1 GHz) and the 
7-ray one (e.g. Fossati et al. 1998, Salamon & Stecker 1996) of 
blazars detected by EGRET. Bloom (2008) found L.y oc L°"±""=' 
with a sample of 122 sources identified as blazars in the revised 
EGRET sample. Mucke et al. (1997) argued that the correlation 
could be due to instrumental biases coupled to the use of average 
7-ray fluxes that washes out the considerable variability of blazars 
at 7-ray wavelengths. Nontheless, Zhang et al. (2001), through par- 
tial correlation analysis, showed that a marginal correlation exists 
between the radio and the 7-ray luminosity in EGRET blazars con- 
sidering the high and low state fluxes of EGRET sources nearly si- 
multaneously observed in the radio band. More recently, Pushkarev 
et al. (2010) showed that there exists a strong correlation between 
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Sources 


L/~y—L/r 


L-,-z 


L, — z 


L.,-Lr (z) 


All (144) 


0.8 
10-34 


0.89 

10-45 


0.72 
10-25 


0.5 
10-11 


FSRQs(112) 


0.66 

10-15 


0.83 

10-30 


0.54 
10-10 


0.44 
10-7 


BL Lacs (22) 


0.78 


0.90 


0.78 


0.29 




10-5 


10-9 


10-5 


10-1 


AGNs (10) 


0.96 


0.97 


0.94 


0.56 




10-^ 


10-^ 


10-5 


10-1 


FSRQs+BL Lacs (134) 


0.76 
10-2^ 


0.86 

io-« 


0.67 
10-18 


0.5 
10-10 



Table 1. Partial correlation analysis of the Lj-Lr correlation accounting 
for the common redshift dependence. Each line gives the Spearman corre- 
lation coefficients and in the last column it is reported the partial correlation 
coefficient. The probabilities of the correlation coefficient are also given. 
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Figure 7. K-corrected 7— ray luminosity (integrated above 100 MeV) ver- 
sus radio luminosity (uLi, at 20 GHz) for the 112 FSRQs, 22 BL Lacs and 
10 "AGNs" (different symbols as shown in the legend) of the AT20G-1LAC 
association sample with measured redshifts. The solid line represents equal- 
ity. The dashed line is the fit considering all the sources (slope 1.13) and the 
dot-dashed line is the fit considering only FSRQs and BL Lacs (slope 1.07). 



the 7-ray and the VLBA radio flux on monthly timescales and that 
the radio flux lags the 7-ray one by 1-8 months. 

Among the 230 AT20G-1LAC associations there are 144 
sources (112 FSRQs, 22 BL Lacs and 10 "AGNs") with measured 
redshifts. In Fig. [7] we show them in the 7-ray versus radio lumi- 
nosity plane. Both luminosities have been k-corrected using the ra- 
dio and 7-ray spectral index of individual sources. While the 7-ray 
luminosity is integrated above 100 MeV, the radio one is the vLi, 
luminosity computed at 20 GHz. FSRQs, BL Lacs and "AGNs" are 
distributed along a linear correlation. Considering only FRSQs and 
BL Lacs (circles and squares in Fig.|7} the correlation has a slope 
1.07±0.05 (dot-dashed line in Fig.|7}. A somewhat steeper slope 



1.13±0.04 is found if also "AGNs" are included in the fit (i.e. con- 
sidering all 144 sourcesjj. This value is steeper than that found by 
Bloom (2008) and also of the correlations (both for the low/high 
states and for the average flux case of EGRET blazars) reported by 
Zhang et al. (2001). 

Given the common dependence of the 7-ray and radio lumi- 
nosity on the redshift z we should test if the L-,-Li correlation 
is true. Several methods have been applied to investigate this pos- 
sibility (e.g. Mucke et al. 1997, Zhang et al. 2001, Bloom et al. 
2008). We perform a partial correlation analysis by removing the 
dependence of both L^ and Li on the redshift z (e.g. Padovani 
1992). We computed the Spearman correlation coefficients and the 
associated probabilities and then the partial correlation coefficient 
and the probability of the null hypothesis that the two luminosities 
are uncorrelated. The chance probability of the partial correlation 
coefficient is distributed as a t-statistic. All the values of the cor- 
relation coefficients and the associated probabilities are reported in 
Tab. 1. We note that considering the FSRQs and BL Lacs together 
the partial correlation probability of the null hypothesis is lO-i". 
This result indicates, in agreement with that reported by Bloom 
(2008), that indeed a strong L^-Li exists in blazars. However, by 
considering BL Lacs and AGNs separately, we find a high chance 
probability of the partial correlation coefficient. This suggests that, 
although we still have few sources of these classes, their large red- 
shift spread makes the luminosity correlation less statistically sig- 
nificant than for the class of FSRQs. 



5.5 Contribution of blazars to the EGBR 

The existence of a L-^-Li correlation and of a corresponding corre- 
lation in the observer frame (i.e. the F^-Fi correlation found in this 
paper) could also have some implications for the computation of the 
contribution of blazars to the extragalactic 7-ray background radi- 
ation (EGBR). One method often adopted to this aim uses a linear 
relation between the 7-ray luminosity of blazars and the luminos- 
ity at some other wavelength in order to re-scale the 7-ray lumi- 
nosity function through the often better known luminosity function 
at the other wavelength (e.g. Salamon & Stecker 1996, Norumoto 
& Totani 2006). Alternatively one can construct the 7-ray lumi- 
nosity function of blazars starting from a catalog, like the Fermi 
first blazar catalog. The latter method has been recently applied by 
Abdoetal. (2010c). 

Fermi finds (Abdo et al. 2010d) that the EGBR spectrum is 
consistent with a power law with spectral index 2.41 ±0.05 and 
an integrated (>100 MeV) flux of l.OSxlQ-'^ cm-^ s-i sr-i, 
which is softer and less intense with respect to the measurements of 
EGRET. Note that the above value is not the total background, but 
the one obtained subtracting out the detected sources. It contains 
the contribution of undetected sources which have either a flux be- 
low the Fermi sensitivity threshold (corresponding to the flux of 
the faintest source detected by Fermi) or that are not detected (but 
with a flux larger than this limit) because of their intrinsic proper- 
ties (e.g. soft spectrum or position in the sky coincident with re- 
gions of high diffuse background level - see Abdo et al. 2010c). 
Recently, Abdo et al. (2010c) considered the contribution of blazars 
to the EGBR. They point out that, due to the detection efficiency of 
Fermi (that we also used in this work), there is a substantial frac- 
tion of 7-ray sources which are not detected but still have a flux 



All the fits are performed with the bisector method (e.g. Isobe et al. 1990). 
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larger than the flux of the faintest source detected by Fermi. There- 
fore, from their count distribution, they estimate that non-detected 
blazars, but with a flux larger than the faintest detected source limit, 
should contribute the ~16% of the EGBR flux (to which detected 
Fermi sources have been subtracted). By extrapolating the blazars' 
count distribution to zero flux, this estimate becomes 23% (Abdo 
et al. 2010c). 

We can perform a simple exercise: we assume that all the 
radio sources in the AT20G survey with flat radio spectrum are 
candidate blazars emitting in the Fermi energy band. Through the 
F-y-Fr correlation we can compute their integrated F^ that can be 
compared with the level of the 7-ray background. This estimate 
should be compared with the EGBR flux including the detected 
sources, i.e. roughly a factor 1.3 larger than the EGBR flux used 
in Abdo et al. (2010c) from which the detected sources were re- 
moved. The diffuse EGBR we use is adapted from Fig. 3 of Abdo 
et al. (2010d) where the EGBR and the contribution of detected 
sources are shown separately. 

We use the reconstructed F^-F-r correlation with slope in the 
range 1.25-1.4 (corresponding to the acceptable solutions in Fig.|6] 
- filled circles) and perform a set of 300 simulations first assigning 
to the 3686 flat radio spectrum AT20G sources an F^ through this 
assumed correlation and then assigning an F-y according to a log- 
normal 7-ray flux variability with o^ — 0.5 (§4). For each simu- 
lation we calculate the total 7-ray flux contributed by the flat radio 
sources. On average we find that their contribution (according to the 
assumed correlation slope and normalization) ranges between the 
~37±5% and the ~52±5% (for slopes 1.25 and 1.4, respectively) 
of the EGBR including the detected sources. This fraction is made 
by two contributions: the total flux of detected sources (i.e. the 230 
AT20G-Ferm( associations) which can be between 20% and 35% 
of the EGBR, while the remaining 17% is the contribution of unde- 
tected sources belonging to the population of flat radio sources with 
Fr^40 mjy (i.e. in the AT20G survey). Based on model population 
studies of blazars, Inoue et al. (2010) finds that the contribution of 
blazars to the EGBR (including detected sources) should be 45% 
(an additional 35% should be due to non-blazar AGNs). This es- 
timate is consistent with the range derived from our analysis. We 
note that in our estimate we are considering the combined contri- 
bution of FRSQs and BL Lacs, although they have different 7-ray 
spectral properties (the latter, having a harder spectrum in the Fermi 
band, are very likely dominating the contribution to the EGBR at 
high energies) and redshift distributions. However, we have used 
the Fry-Fr correlation for our estimate and we have shown (GIO) 
that both these classes of sources follow a similar correlation be- 
tween the radio and the 7-ray flux. 

5.6 Predictions for the brightest 7-ray blazars 

One possible application of the F^-F^ correlation obtained in this 
paper is to predict the average 7-ray flux and its maximum value 
for any given blazar with known radio flux density at 20 GHz. Fig.[8] 
shows the AT20G-1LAC associations (GIO) and the reconstructed 
correlation (here we have chosen to report the solution with slope 
1.25) and its 1,2 and 3 a^ dispersion. The F^-F^ correlation di- 
vides the plane of Fig.[7]into two regions that we label as "low" and 
"high" state. These correspond to those states of the 7-ray flux (1,2 
or 3 a) above or below the average value represented by the F-y-Fi 
correlation (solid line in Fig.[8}. The shaded regions in Fig.[8]show 
how the average flux (over 1 1 months as measured by Fermi in its 
survey which is the one we used to construct the F^-Fr correlation) 
can vary according to a log-normal distribution with a^ ~0.6. 



We have considered the brightest blazars, with Fr at ~ 20 
GHz larger than 3 Jy, distributed in the southern (i.e. present in the 
AT20G survey) and northern emisphere. For them we can calculate 
the average _P-, and the maximum average F-y they can reach if their 
long-term variability follows the log-normal distribution found for 
the sources in common between EGRET and Fermi. These values 
of F.y and F.y,scr are given in Tab. |2] Among the sources in which 
the 7-ray flux can be larger than 5 x 10~^ ph cm~^ s~^ there are 
3C 279 and 3C 273. Furthermore, Tavecchio et al. (2010) showed 
that in the brightest blazars a short-term variability (§3) can be 
present and modulate the 7-ray flux by a factor 3 (or even 10 in 
the most extreme cases), on top of the long term variability. There- 
fore, they could reach even larger fluxes through sporadic flares of 
emission, like the case of 3C 454.3 last year (Bonnoli et al. 2010; 
Foschini et al. 2010), increasing their Ky from the highest average 
value (i.e. the F.y^3cr) still by a factor of a few. Tab. |2] provides the 
list of those sources which can be included in long-term multi- 
wavelength monitoring of blazars for the study of their variability 
and for the characterization of the most extreme phases of their 
emission. 

An interesting case is represented by 4C-I-21.35 (PKS 
1222-1-216 at z=0.43). In April 2009 this source increased its av- 
erage Fj flux by a factor of 10 with respect to its average flux of 
~4.6x 10~* phot cm~^ s~^ measured in the Fermi first 6 months 
survey (Longo et al., 2009). A strong flare was detected by AG- 
ILE/GRID in December 2009 (Verrecchia et al. 2009) with a flux 
(integrated above 100 MeV) of 2.5 x 10~® phot cm"^ s^^ and con- 
firmed by Fermi (with flux 3.4x10"^ phot cm~^ s~^ - Ciprini et 
al., 2009). In the period April-May 2010 GeV flares were detected 
by Fermi (Donato et al., 2010) at a flux level of 8x 10"'' phot cm~^ 
s~^, i.e. a factor about 4 in excess with respect to the average flux 
of that period and by AGILE (Bulgarelli et al., 2010). Very High 
Emission (above 100 GeV) was also found from this source in this 
period (Neronov et al., 2010; Mariotti et al., 2010). Finally in June 
2010 Fermi recorded a flux of 1.2x 10^'' phot cm~^ s^^ (lafrate et 
al., 2010) which represents an increase of a factor of 3 with respect 
to the average flux of the week. 

4C-I-21.35 has a radio flux (at 15 GHz) of about 1 Jy. The F^~ 
Fr correlation can be used to calculate its average flux level and 
its average 3a" flux level. These turn out to be F^=2.2x 10~* phot 
cm~^ s~^ and ^^^3^=1.4x10"^ phot cm~^ s"-'. We note that the 
brightest flare of 4C-I-21.35 detected by Fermi in June 2010, has a 
flux a factor ~10 larger than the maximum average flux predicted 
from the F^-F^ correlation, i.e. if the source were 3a brighter than 
its average flux as predicted by the correlation. 

As already discussed, one possibility is that while the average 
flux (i.e. the average as measured by Fermi in its 1 1 months survey) 
is modulated by a long-term variability following a log-normal 
distribution as found in the case of the EGRET sources detected 
by Fermi (see §3), on top of this there is a much shorter variability 
(described by Eq.lTJ as found in the brightest blazars (Tavecchio et 
al. 2010). This short-timescale variability operates on top of each 
state of the average 7-ray flux and can boost the flux still by a fac- 
tor 3-10 even when the source has reached its maximum value of 
the average flux. 

In Fig. [8] we show the variations of the 7-ray flux of PKS 
1222-1-216 caught to be in outburst by both Fermi and AGILE on 
December 2009 (Verrecchia et al. 2009, Ciprini et al. 2009) and 
to reach in June 2010 the flux of 1.2x10"^ phot cm~^ s~\ We 
also report the recent detected outburst of PKS 1830-21 (Ciprini et 
al. 2010) which went through an outburst in October 2010 reach- 
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Figure 8. 7-ray flux vereus radio flux density of the AT20G-1LAC associa- 
tions (open red squares). The solid line represents the reconstructed F-y-Fr 
correlation and its 1,2 and 3 cr-y scatter (with a^y = 0.6) is shown by the 
shaded regions. The names of the brightest sources are shown. The region 
above the solid line (with slope 1.25) represents the high state: i.e. a source 
can vary its average F^ within the shaded regions according to a log-normal 
probability function with assigned cr. The dashed line is the Fj-Fr corre- 
lation. The outbursts detected by Fermi of two sources, i.e. PKS 1222-1-216 
and PKS 1830-21, are shown. 



ing a daily F^ of 5.2x10 ^ phot cm ■^ s ^ with a peak flux of 



1.4x10" photcm"^ s" 



6 DISCUSSION AND CONCLUSIONS 

We studied the F^-Fr correlation between the 7-ray (> 100 MeV) 
flux and the radio flux (at 20 GHz) observed in the 230 Fermi 
sources with a counterpart in the 20 GHz ATCA survey (GIO). This 
correlation is biased by the radio flux limit of the AT20G survey and 
at low 7-ray fluxes by the Fermi sensitivity. However, in the F-y- 
Fr plane there are regions (hatched triangles in Fig.[TJ where these 
selection effects are not present and still no source is found. This 
suggests that the F-/~Fr correlation is real. And yet, only 1/15 of 
the radio sources with flat radio spectra in the AT20G survey have 
a counterpart in the 1 1 months Fermi survey. 

Through numerical simulations we have recovered the true 
F-y-F,- correlation that can reproduce the observed F^-F,- one. In 
doing this we have considered the two main instrumental selection 
effects (radio flux limit and Fermi sensitivity) and we tested the 
possibility that the non detection of radio sources by Fermi could 
be due to the variability of their emission in the 7-ray energy range. 
Tavecchio et al. (2010) characterized the variability of 3C 454.3 and 
PKS 1510-089, two among the brightest blazars detected by Fermi 
(see also Foschini et al. 2010). They found a daily flux distribu- 
tion which we here modeled through Eq. [T] However, this "short" 
timescale variability induces a small variation of the 7-ray flux 
when we average it over 1 year. A longer timescale variability is 
found by comparing, for the sources in common, the fluxes mea- 
sured by EGRET and, after almost ten years, by Fermi. The ratio of 
the 7-ray fluxes of the blazars detected by EGRET and Fermi have 



Source 


Fr 


F., 


P'j/iT 




■fy 


phot cm ^ s ^ 


phot cm^^s^^ 


3C279 


20.0 


9.5e-7 


6.0e-5 


[HB89] 1921-293 


13.8 


6.0e-7 


3.8e-5 


3C446 


8.3 


3.2e-7 


2.0e-5 


PKS 0727-11 


6.7 


2.4e-7 


1.5e-5 


PKS 1335-127 


6.1 


2.1e-7 


1.4e-5 


PKS 0420-01 


6.0 


2.1e-7 


1.3e-5 


PKS 1830-21 


5.5 


1.9e-7 


1.2e-5 


PKS 0537-441 


5.3 


1.8e-7 


l.le-5 


PKS 0454-46 


4.2 


1.3e-7 


8.5e-6 


CRATES J0609-0615 


4.2 


1.3e-7 


8.5e-6 


PKS 0402-362 


4.0 


1.2e-7 


8.0e-6 


PKS B0607-157 


4.0 


1.2e-7 


8.0e-6 


PKS 0521-36 


4.0 


1.2e-7 


8.0e-6 


PKS 0454-234 


3.8 


1.2e-7 


7.5e-6 


PKS 1954-388 


3.8 


1.2e-7 


7.5e-6 


APLIB 


3.4 


l.Oe-7 


6.5e-6 


PKS 0208-512 


3.3 


l.Oe-7 


6.3e-6 


PKS 2227-08 


3.2 


9.6e-8 


6.0e-6 


PKS 0637-75 


3.1 


9.2e-8 


5.8e-6 


PKS 1510-08 


2.9 


8.5e-8 


5.3e-6 


3C273 


23.8 


1.2e-6 


7.4e-5 


3C345 


12.0 


5.0e-7 


3.2e-5 


3C 454.3 


11.0 


4.5e-7 


2.3e-5 


87GB[BWE91] 0059-H5808 


8.6 


3.3e-7 


2.1e-5 


[HB89] 2145-1-067 


8.5 


3.3e-7 


2.1e-5 


OJ -1-287 


6.0 


2.1e-7 


1.3e-5 


[HB89] 0735-1-178 


5.3 


1.8e-7 


l.le-5 


[hB89] 2134-1-004 


5.1 


1.7e-7 


l.le-5 


iHB89] 0923-1-392 


5.0 


1.7e-7 


l.le-5 


BLLac 


4.5 


1.5e-7 


9.2e-6 


[HB89] 2201-1-315 


4.5 


1.5e-7 


9.2e-6 


4C -1-50. 11 


4.2 


1.3e-7 


8.4e-6 


[HB89] 1055-1-018 


4.2 


1.3e-7 


8.5e-6 


[hB89] 1308-)-326 


3.9 


1.2e-7 


7.6e-6 


[hB89] 1611+343 


3.6 


l.le-7 


6.9e-6 


[hB89] 1928+738 


3.5 


l.le-7 


6.7e-6 


LBQS 0106+0119 


3.5 


l.le-7 


6.8e-6 


[HB89] 2005+403 


3.4 


l.Oe-7 


6.4e-6 


[hB89] 0234+285 


3.4 


l.Oe-7 


6.5e-6 


[hB89] 0642+449 


3.3 


l.Oe-7 


6.4e-6 


iHB89] 1749+096 


3.3 


l.Oe-7 


6.3e-6 



Table 2. Blazars with the highest radio flux density at 20 GHz. i.e. Fr^3 
Jy, distributed in the southern and northern emisphere (top and bottom table, 
respectively). For each source it is reported its name, the radio flux density 
(for southern sources this is extracted from the AT20G survey - Murphy et 
al. 2009) while for the northern sources it is taken from the Nasa Extragalac- 
tic Database and in most cases it is at 22 GHz. We give the average 7-ray 
flux i^-, calculated from the reconstructed F-y-Fr congelation (i.e. solid line 
in Fig. [S) and its 3 O" highest value con'esponding to the highest possible 
state according to a log-normal variability of the 7-ray flux with cr-y=0.6 
(shaded region in Fig.[8). 



a log-normal distribution with a standard deviation a-, = 0.5, and 
12 of the 66 blazars detected by EGRET have not been detected by 
Fermi. This suggests that the "short" timescale variability is super- 
imposed on a "decadal" variability and that the latter has a larger 
relative variation. We considered these two variability patterns in 
our simulations aimed at recovering the true F-,-Fv correlation of 
the AT20G-Fermi associations. 

As shown in Fig. [5] we cannot reproduce the observed F^~Fr 
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correlation if we adopt the short timescale variability (and average 
the flux over 1 year). In this case we over predict the number of 
simulated sources detectable by Fermi for any combination of the 
F-y-F-c correlation slope and normalization. Instead, we find (Fig. 
|6ll that the distribution of real sources in the F-,-Fi plane can be 
reproduced if the slope of the true F-,-Fr correlation is in the range 
1. 3- 1. 5 and if the 7-ray (1-year averaged) flux variability is mod- 
eled as a log-normal with a^ Js 0.5 (acceptable solutions in Fig.|6] 
represented by filled and open circles). We also find that a smaller 
flux variability, also modeled as a log-normal, cannot reproduce the 
observed F-,-Fi correlation, similarly to what found by assuming 
the short timescale variability pattern of Eq.[T] 

Therefore, we need a "decadal" 7-ray variability which mod- 
ulates the I-year averaged 7-ray flux by at least a factor 3 in order 
to explain the radio/7-ray statistical behavior of blazars. The need 
to assume such a variability pattern is corroborated by the findings 
on the ratio of EGRET to Fermi flux for the sources in common, 
although this is based on the EGRET sources which are the most 
luminous (and maybe the most variable). However, our results show 
that at least a variability of a factor 3 in flux is necessary to repro- 
duce the observed radio-7 flux correlation. Recent results, obtained 
after the 1 1 months Fermi survey, are illuminating: PKS 1222-1-216 
increased its average 7-ray flux by a factor 300 with respect to its 
average flux during the first 6 months Fermi survey. This implies 
that not only the 7-ray flux can vary on very short timescales, but 
also on longer timescales and with an amplitude that is larger than 
what seen in the radio. 

For our simulations, we we used the Fermi sensitivity com- 
puted by Abdo et al. (2010c) which assumes a spectral index dis- 
tribution typical of FSRQ sources. Although the Fermi sensitiv- 
ity depends on the source spectral index (e.g. Abdo et al. 2010a), 
our sample of 1FGL-AT20G associations is dominated by FSRQ 
sources. Nonetheless, for BL Lac objects, due to their harder spec- 
trum, a lower Fermi detection limit (represented by the shaded re- 
gion S2 in Fig[T]l should be considered for these sources. We re- 
mark that based only on the radio data we cannot distinguish in 
our simulations between FSRQ and BL Lac objects among the ra- 
dio sources with flat radio spectrum but without a 7-ray counter- 
part. However, if a lower detection sensitivity would be used in 
our simulations, a combination of a slightly smaller normalization 
and slope of the F-y-F^ correlation possibly coupled with a slightly 
larger variability of the 7-ray flux would be required to reproduce 
the F-y-Fi correlation. Still the results would be comparable with 
those derived with the detection sensitivity of FSRQ. 

As a result of the existence of the F^-Fr correlation and of the 
long-timescale variability we also expect that the brightest radio 
blazars (e.g. with radio flux density Js3 Jy) will always be detected 
by Fermi since the variability of the 7-ray flux cannot make them 
much fainter (on average, short duration events at very low fluxes 
can occur). Our result also implies that if the average fluxes in the 
radio and 7-ray bands are correlated, then also a decadal radio flux 
variation should be expected and with a positive time delay if the 
radio emitting region is larger than the one emitting the 7-rays, and 
thus is located further out along the jet. 

However, our fluxes are collected at different epochs: the radio 
data are single snapshots obtained by the AT20G survey in the pe- 
riod 2004-2008 while the Fermi fluxes are averages over 1 1 months 
survey between 2008 and 2009. Pushkarev et al. (2010) have shown 
that there exists a strong correlation between the 15.4 GHz core 
VLBA flux density and the 7~ray one and that these two fluxes are 
correlated with a typical delay distributed between 1 and 8 months 



(in the observer frame) which is interpreted as caused by opacity 
effects. 

This study, on the correlation of fluxes, is related to the possi- 
bility that there exists a correlation between the 7-ray and the ra- 
dio luminosities in blazars (as found e.g. by Bloom et al. 2008 and 
discussed in Mucke et al. 1997). Although based on the fraction 
(144) of the AT20G-Ferm; associations with measured redshifts, 
we showed that a strong correlation exists between the radio lu- 
minosity and the 7-ray one and that, considering FSRQs and BL 
Lac objects, it is linear. We also verified, through partial correlation 
analysis, that this correlation is not due to the common dependence 
of the luminosities on redshift for FSRQs (null hypothesis proba- 
bility of the partial correlation, removing the redshift dependence, 
P = 10~^"), while only a marginal claim can be made for BL Lacs 
alone. 

The other main consequence of the existence of a F-y-Fr cor- 
relation in blazars is that it could be used to estimate the contribu- 
tion of these sources to the 7-ray background (Abdo et al. 2010c). 
The contribution of blazars can be estimated considering the true 
Fry-Fr correlation and the population of radio sources with flat ra- 
dio spectra (i.e. candidate blazars) in the AT20G survey. We esti- 
mate that the blazar contribution to the extragalactic diffuse back- 
ground is between 37% and 52% (according to the parameters of 
the assumed F-y-Fr correlation) of which 17% is the contribu- 
tion from non-detected sources. Considering that our estimates are 
based on a flux limited sample of radio sources (i.e. those with 
20 GHz flux larger than 40 mJy) we should expect that if the ra- 
dio flux limit is further decreased, the number of sources should 
increase. While detailed predictions depends on the still poorly 
known LogN — LogS of radio sources at very low flux levels, 
we should expect that the contribution of blazars to the EGRB can 
even be larger than our present estimates based on a radio flux lim- 
ited sample. 
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